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Glycinergic inhibition tunes coincidence detection 
in the auditory brainstem 
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Neurons in the medial superior olive (MSO) detect microsecond differences in the arrival 
time of sounds between the ears (interaural time differences or ITDs), a crucial binaural cue 
for sound localization. Synaptic inhibition has been implicated in tuning ITD sensitivity, but 
the cellular mechanisms underlying its influence on coincidence detection are debated. Here 
we determine the impact of inhibition on coincidence detection in adult Mongolian gerbil 
MSO brain slices by testing precise temporal integration of measured synaptic responses 
using conductance-clamp. We find that inhibition dynamically shifts the peak timing of 
excitation, depending on its relative arrival time, which in turn modulates the timing of best 
coincidence detection. Inhibitory control of coincidence detection timing is consistent with 
the diversity of ITD functions observed in vivo and is robust under physiologically relevant 
conditions. Our results provide strong evidence that temporal interactions between excitation 
and inhibition on microsecond timescales are critical for binaural processing. 
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Coincidence detection of afferent synaptic inputs is a 
hallmark of neuronal processing. In many brain areas, 
such as the cerebral cortex and cerebellum, neurons 
integrate hundreds of synaptic inputs and perform coincidence 
detection with millisecond resolution. In contrast, neurons in the 
mammalian brainstem nucleus of the medial superior olive 
(MSO) integrate relatively few, clearly defined synaptic inputs^ 
but must perform coincidence detection with orders of magnitude 
higher precision. They detect microsecond differences in the 
arrival time of low- frequency sounds between the ears (interaural 
time differences or ITDs)^ and generate a rate code that 
represents the location of sounds in the horizontal plane^. 

MSO neurons generally respond best to contralateral-leading 
ITDs"*"^. To allow best coincidence detection of such temporally 
staggered inputs, the circuit must generate internal delays to 
compensate for the external delay created by the head. A longer 
excitatory path length from the contralateral side can explain this 
phenomenon at the population level, but there is increasing 
evidence that axonal length disparity^ alone is insufficient to 
explain the diversity of preferred ITDs observed in vivo^^~^^. 
Thus, it has become a crucial question how individual neurons 
within the population are tuned to their specific ITDs^'^"*"^^, an 
answer to which is necessary for a complete understanding of 
low-frequency sound localization in mammals. 

Previous studies found that blocking glycinergic inhibition 
in vivo shifted the peak of ITD functions, providing direct 
evidence that inhibition is important for ITD tuning^'^. 
As glycinergic inputs to the MSO are extremely fast and 
precise^ , a computational model was developed in which 
the relative arrival time of excitatory and inhibitory inputs tunes 
neurons to their preferred ITDs by modulating the peak timing of 
the summated input from each side^'^^'^"*'^^. However, this model 
was not rigorously tested under naturally occurring conditions 
and has been recently challenged^'^^. In one study, the difference 
in arrival time between inhibition and excitation was measured 
in vitro. Inhibition was indeed precise and actually preceded 
excitation, but the average timing condition did not tune ITD 
functions, and it was concluded that inhibition was not important 
for ITD tuning^^. In a second study, ITD functions were 
determined by juxtacellular and whole- cell recordings in vivo, 
but preceding inhibition as measured in vitro could not be 
detected^. Thus, there are now discrepancies between the existing 
data^"^'^"*"^^ and confusion regarding the relevance of inhibition 
for ITD tuning^'2^-29. 

To understand these discrepancies, we systematically re- 
evaluate the role of precise inhibition in synaptic integration 
and coincidence detection under well-defined, physiologically 
relevant conditions using acute brain slices from adult Mongolian 
gerbils. We find that inhibition indeed tunes coincidence 
detection in a manner consistent with the original model^'^'^"*'^^, 
but much more dynamically than originally proposed. Moreover, 
we provide evidence that the specific involvement of inhibition in 
ITD tuning in vivo may not be discernable with currently 
available techniques. Thus, we expand the framework for the role 
of inhibition in ITD tuning and provide explanations for the 
discrepancies in the literature. 

Results 

Inhibition modulates the peak timing of excitation. Coin- 
cidence detection neurons in the MSO are strikingly bipolar in 
shape (Fig. la)^^ and receive bilateral excitatory and inhibitory 
synaptic inputs (Fig. Ib)^^ Because synaptic kinetics and input 
resistance change markedly during development^^'^^'^^, we 
restricted our primary experiments exclusively to mature gerbils 
of ages postnatal day (P) 60-90. As a basis for our quantitative 



re- evaluation of cellular ITD sensitivity, we first extracted the 
time course of excitatory and inhibitory synaptic conductances 
(EPSGs and IPSGs, respectively) in a separate set of 
voltage-clamp experiments (Supplementary Fig. 1). Individual 
conductance waveforms were then selected and injected as 
templates into the soma using conductance -clamp to simulate 
excitatory and inhibitory postsynaptic potentials (EPSPs 
and IPSPs, respectively), which produced events that largely 
resembled synaptically evoked responses (Supplementary Fig. 2). 
This allowed us to investigate precise temporal interactions 
between excitation and inhibition based on synaptic responses 
measured in mature tissue. 

A fundamental principle of the original inhibitory model of 
ITD tuning proposes that inhibition modulates the timing of best 
coincidence by enforcing a shift in the peak timing of excitation 
on each side^'^'^"^'^^. As there is strong evidence that the arrival 
times of excitatory and inhibitory inputs to the MSO is exquisitely 
precise^^'^^'^"*, we determined whether a precisely timed IPSP can 
shift the peak timing of an EPSP, here arbitrarily designated 
to the contralateral side (Fig. 1). Inhibitory inputs are 
predominantly restricted to the soma^^~^^, and thus simulating 
IPSPs with conductance- clamp is suitable to mimic realistic 
IPSPs. However, because excitatory inputs predominantly 
target the dendrites^^ we compared the influence of a single 
conductance -clamp -simulated (G-inject) IPSP on the peak 
timing of a synaptically evoked (F-stim) EPSP (Fig. Ic) and a 
conductance -clamp -simulated EPSP (Fig. Id) in the same 
recording. We used an IPSG template with a decay time 
constant (tipsg) that represented our population average 
(1.5 ms) and an EPSG template with a decay time constant 
('^^epsg) that generated an EPSP matching the kinetic profile of the 
synaptically evoked EPSP (0.2-0.3 ms, where 0.3 ms was the 
average; Supplementary Fig. 2). EPSG and IPSG peak amplitudes 
were then adjusted to achieve an effective excitation-inhibition 
(E-I) ratio of 1:1 (see Methods). 

Because the actual arrival times of excitation and inhibition to 
the MSO in vivo is not known, we investigated a broad range of 
relative inhibitory timing conditions (Atinh; Fig. le). Indeed we 
found timing conditions that advanced the peak timing of both 
the synaptically evoked and conductance- clamp -simulated EPSPs 
(Atinh = 0.1 ms, F-stim: — 62 ± 8 jis; G-clamp: — 53 ± 3 jis, n — 7\ 
Fig. lf,g, left) and others that delayed the peak of the EPSPs, albeit 
to a lesser extent (Atinh= —0.6 ms, F-stim: 21 ±9 jis; G-clamp: 
4±6|is, n = 7; Fig. lf,g, right). The similarity of the timing- 
dependent peak shifts between synaptically evoked and con- 
ductance-clamp-simulated EPSPs (Atinh = 0.1 ms: P = 0.986; 
^^inh = — 0-6 ms: P — 0.647 , two-way analysis of variance 
(ANOVA)) indicates that inhibition influences synaptic integra- 
tion at the soma (Fig. Ih). This was further supported by separate 
dual- electrode conductance-clamp experiments that revealed no 
substantial difference in PSP kinetics or the magnitude of 
inhibition- enforced EPSP peak shifts compared with single- 
electrode experiments (Supplementary Fig. 3). Taken together, 
these findings demonstrate that single- electrode conductance- 
clamp is suitable to investigate precise temporal interactions 
between excitation and inhibition in MSO neurons. Thus, 
inhibition can indeed modulate the peak timing of excitation 
from either side in a bidirectional manner depending on its 
arrival time relative to excitation. These findings are consistent 
with the basic predictions of the original inhibitory model^'^'^"*'^^ 
as well as recent reports^^. 

Input speed and balance influence peak shift magnitude. 

Conductance-clamp is advantageous over synaptic stimulation in 
that many parameters can be tested quickly in a single recording. 
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Figure 1 | Inhibition modulates the timing of excitation, (a) Fluorescence micrograph of a P60 MSO neuron. Scale bar, 50|im. 
(b) Schematic of input circuitry to the MSO. Ipsilateral (Ipsi) and contralateral (Contra) glutamatergic inputs from the anteroventral cochlear nuclei 
(AVCN) target the dendrites. Glycinergic inputs from the lateral and medial nucleus of the trapezoid body (LNTB and MNTB, respectively) target the 
soma. (c,d) Schematics of configurations for synaptic stimulation of EPSPs (F-stim) while simulating IPSPs (G-inject) with conductance-clamp (G-clamp, c) 
and simulating both EPSPs and IPSPs with conductance-clamp (d). (e) Timing convention. The onset of Inhibition is timed relative to that of Excitation from 
one side (Atinh), designated as contralateral. Negative values indicate Inhibition leads Excitation (Inh-lead), and positive values indicate Inhibition lags 
Excitation (Inh-lag). (f,g) Voltage traces for an example recording using configurations in c,d, respectively, aligned in time to the peak of the EPSP alone 
(light traces). The influence of a precisely timed IPSP (red traces) on EPSP peak timing is visualized as a shifted peak of the composite PSP (dark traces). 
Traces are separated to illustrate inhibitory timing conditions that enforced a peak advance (left) or delay (right). Traces are also plotted aligned in time to 
the IPSP in Supplementary Fig. 2e. Resting membrane potential (Kest): — 64mV. (h) Average (±s.e.m.) EPSP peak shifts plotted against Atjnh for 
synaptically evoked (blue) and conductance-clamp-simulated (black) EPSPs. Negative and positive values indicate an EPSP peak advance and delay, 
respectively. Atjnh = 0.1 ms, P = 0.986; Atjnh = -0.6 ms, P = 0.647; two-way ANOVA, n = 7 recordings. 



Using this advantage, we investigated a broad parameter space in 
each neuron to gain insight into the mechanism by which the 
relative timing of inhibition modulates the peak timing of 
excitation. For example, this allowed us to test whether the wide 
range of excitatory and inhibitory input kinetics we observed 
(Supplementary Fig. 1) is relevant for the efficacy of inhibition- 
enforced EPSP peak shifts (Fig. 2a-f). Here, we compared 
representative fast (Fig. 2a,b) and slow (Fig. 2c,d) EPSG and IPSG 
templates with their average-speed (tepsg = 0-3 ms, Tipsg=1-5 
ms) counterparts (Supplementary Fig. 2), maintaining an E-I 
ratio of 1:1. Peak shifts were substantially dependent on excitatory 
input kinetics, where slower excitation (tepsg = 0-5 ms) was 
comparably more sensitive to peak shifts in both directions 
(Atinh = 0.1 ms: - 92 ± 12 versus - 58 ± 6 |is, P = 0.003; Ati^h = 
-0.6ms: 32 ± 8 versus 2±6)is, P = 0.011, two-way ANOVA, 
n = 5; Fig. 2a,c,e). Peak shifts were also sensitive to inhibitory 
input kinetics, where faster inhibition (tipsg= 1.0 ms) was more 
effective at enforcing the peak delay (Atinh = — 0.6 ms: 20 ± 6, 
P = 0.015, two-way ANOVA, n = 5), but did not substantially 
impact the peak advance (Atinh = 0.1 ms: — 62 ± 8, P = 0.925, 
two-way ANOVA, n = 5; Fig. 2b,d,f). Thus, although we 
predominantly use the average-speed EPSG and IPSG throughout 
the manuscript, slower excitation and faster inhibition will 
generally enhance the impact of inhibition on excitatory timing. 

Because inhibition can impose both a shunt and hyperpolar- 
ization on the membrane, we next evaluated their independent 
contributions to EPSP peak shifts (Fig. 2g). Shunting inhibition 



was isolated by hyperpolarizing neurons to the set reversal 
potential for chloride (-90 mV)^. This manipulation revealed 
a peak advance regardless of timing (Atinh = 0.1 ms: — 32 ± 4 |is; 
Atinh =- 0.6 ms: -23±4|is, n = ll; Fig. 2g, dotted black). 
Depolarizing neurons from this point revealed the contribution of 
hyperpolarization that manifested itself as a biphasic function 
superimposed on the shunt-induced function. Compared with the 
median resting potential ( — 65 mV, Fig. 2g, maroon and 
Supplementary Fig. 4) more depolarized potentials ( — 45 mV, 
black) revealed larger peak shifts in both directions 
(Atinh = 0.1ms: - 101 ± 3 versus -67±4jis, P<0.001; 
Atinh= -0.6ms: 17±6 versus - 10 ± 5 jis, P<0.001, two-way 
ANOVA, n = ll; Fig. 2g). We also doubled the E-I ratio to 
reduce the relative strength of inhibition and found that peak 
shifts were proportionally reduced, particularly for the peak 
advance (at -65mV, Atinh = 0.1ms: -42±4|is, P = 0.002; 
Atinh =- 0.6 ms: -3±6|is, P = 0.972, two-way ANOVA, 
n = ll; Fig. 2h). Thus, an inhibitory shunt generally advances 
the peak of an EPSP, but hyperpolarization enforces the 
bidirectional effect of inhibition on EPSP peak timing, 
depending on the balance of excitation and inhibition. Here, 
more depolarized membrane potentials (imposing a larger driving 
force on IPSPs) and stronger inhibition (relative to excitation) 
enhance the inhibitory control of excitatory timing. 

Because MSO neurons integrate their synaptic inputs 
linearly^'^^, these findings can be explained by a simple model 
(Fig. 2i,j). Here, we schematized an EPSP as a Gaussian function 
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and the hyperpolarizing profile of an IPSP as a sloped line (in this 
case the re-depolarizing phase). If we add the line to the Gaussian 
function, the peak of the summated function is shifted 
proportionally towards the slope of the line (Fig. 2i). This 
explains why coincident or slightly lagging inhibition (imposing a 
negative slope) enforces a peak advance, and why leading 
inhibition (positive slope) enforces a peak delay. It also 
illustrates why faster and stronger inhibition (steeper slope) is 



more effective at enforcing EPSP peak shifts^^. Because 
differences in inhibitory kinetics predominantly influence the 
slope of the re- depolarizing phase of the IPSP (Supplementary 
Fig. 2), it further shows why inhibitory kinetics predominantly 
influenced the peak delay compared to the peak advance. Finally, 
if we double the Gaussian half- width to mimic a relatively slower 
EPSP, larger shifts are generated (Fig. 2j). This explains why 
slower excitation is more sensitive to peak shifts. Thus, inhibitory 
control of excitatory timing is influenced by factors beyond the 
relative time of arrival^'^'^^'^^. 
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The activity of potassium channels facilitates peak shifts. Before 
investigating the consequences of inhibition-enforced EPSP peak 
shifts on bilateral coincidence detection, we investigated the 
role of low-threshold potassium channels (Klta) in inhibition- 
enforced EPSP peak shifts. Klta has been implicated in 
sharpening EPSPs^^'^^ and has recently been shown to interact 
with IPSPs^^. Simply blocking Klta does not allow an assessment 
of its influence on synaptic integration because Klta also 
maintains the neuron's exquisitely fast membrane time constant 
in the adult (~180|is). To circumvent this pharmacological 
limitation, we generated a point-neuron model based on our 
electrophysiological measurements and performed cell-wise 
fittings of EPSPs and IPSPs that were recorded in conductance- 
clamp (Fig. 3). The model recapitulated the kinetic profile and 
amplitude of measured events (Fig. 3a) by altering only a few 
independent parameters (Fig. 3b,c), demonstrating its accuracy 
for evaluating synaptic integration. It also successfully 
recapitulated inhibition -enforced EPSP peak shifts for 
individual recordings (Fig. 3d) and the entire population (Fig. 3e). 

To determine whether native Klta activity influenced 
inhibition- enforced EPSP peak shifts, we eliminated the Klta 



Figure 2 | Input speed and balance influence peak shift magnitude. 

(a-d) Nornnalized (nornn) voltage traces for an example recording of 
inhibition-enforced EPSP peak shifts connparing fast (tepsg = 0.2 ms, 
a; T^IPSG^I-Oi^s, b; goldenrod) and slow (tepsg = 0.5 ms, c; T|psg = 2.2 ms, 
d; chartreuse) speed EPSGs (a,c) and IPSGs (b,d) with their average-speed 
counterparts (tepsg = 0.3 ms; T|psG = 1-5ms) for Atinh = 0.1 (left) and 

- 0.6 ms (right). Insets are zooms of the peaks, aligned in amplitude. Inset 
scale bar (a, bottom), 0.1 mV. Traces are colour-coded as follows: light 
traces indicate the EPSP and IPSP alone, dark traces indicate the composite 
PSP. Traces are aligned in time to the peak of the EPSP alone for each kinetic 
template pair. V^est- -62mV. (e,f) Average (±s.e.m.) EPSP peak shifts 
plotted against Afjnh for comparing EPSG (Afinh = 0.1 ms: tepsg = 0.2 ms, 

P = 0.542; Tepsg = 0.5 ms, P = 0.003; Atjnh = - 0.6 ms: tepsg = 0.2 ms, 
P = 0.422; Tepsg = 0.5 ms, P = 0.011) (e); and IPSG (AU = 0.1 ms: 
T|psG = 1-0ms, P = 0.925; T|psg = 2.2 ms, P = 0.213; Atjnh = -0.6 ms: 
T|psG = 1-0ms, P = 0.015; T|psg = 2.2 ms, P = 0.500) (f) kinetics. Two-way 
ANOVA, n = 5 recordings. (g,h) Same as in e,f, except the average-speed 
EPSG and IPSG were used, and the holding potential (V^) was between 

— 90 and — 45mV in 5mV increments. Compared with the median 
V^est- -65mV (maroon; Supplementary Fig. 4b), -45mV (black): 
Afinh = 0.1ms, P< 0.001; Atj^h^ -0.6 ms, P< 0.001. The E-l ratio was 1:1 
in e, but 2:1 in f. Compared with E-l ratio = 1:1, E-l ratio = 2:1 at 

V^ = -65mV: Atinh = 0.1ms, P< 0.001; Atinh= -0.6 ms, P = 0.972. 
Two-way ANOVA, n = 11 recordings. (i,j) Model illustrating the mechanism 
underlying the dynamics of Inh-enforced EPSP peak shifts. Fast and slow 
EPSPs are modelled as Gaussian functions of comparably narrow (i) and 
wide (j) half-widths, respectively. Slow, weak (pink) and fast, strong (red) 
IPSPs are modelled as shallow and steep sloped linear functions, 
respectively. The EPSP alone is shown in light grey with a peak at t = 0. 
The composite PSP is the linear sum of the Gaussian function and the 
shallow (grey) or steep (black) linear functions. AU, arbitrary units. 
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activation kinetics, thereby rendering the channel passive. 
Interestingly, this manipulation did not substantially reduce peak 
shifts (Atinh = 0.1 ms: — 66 versus — 73 jis; A^inh = — 0.6 ms: — 6 
versus -6 jis; Fig. 3f, turquoise). However, optimal model fit 
parameters (Fig. 3b) supported by measured intrinsic membrane 
properties (Supplementary Fig. 4) indicate that in contrast to 
observations in juvenile (PI 6- 19) gerbils^^ Klta channels in 
adult (P60-P90) gerbils are activated at more hyperpolarized 
potentials (1^0.5= — 73mV), thus rendering a larger fraction of 
them open at rest ( — 65 mV). Indeed, we found that in contrast 
to P17 gerbils, the slowing and tightening of the membrane with 
increasing hyperpolarization (indicating that Klta had closed) 
was not complete, even at -90 mV (Supplementary Fig. 4d,e). 
Interestingly, depolarizing the activation voltage in the model to 
levels observed in juvenile gerbils (Vo.5= — 38mV)^^ revealed 
more robust peak shifts, particularly for a peak delay 
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(Atinh = 0.1ms: - 90 jis; Atinh= -0.6 ms: 10 jis; Fig. 3f, black). 
We also repeated the same protocol using the slow-speed EPSG 
and found even more pronounced peak delays (Atinh = — 0.6 ms: 
55 versus 8 jis; Fig. 3g). 

To gain insight into the mechanism underlying the facilitation 
of Klta activity on inhibition-enforced EPSP peak shifts, we 
further evaluated the EPSPs (Fig. 3h,i) and the IPSP (Fig. 3j) 
produced by the model. When the Klta activation voltage was 
depolarized, both EPSPs and the IPSP amplitudes increased, but 
only IPSP kinetics were visibly altered. Klta gating is fast enough 
to follow the profile of the IPSP and is thus closed during the 
hyperpolarizing phase. This naturally moves the membrane 
potential away from the potassium reversal potential, speeding 
the re- depolarizing phase, from which our previous experiments 
would predict larger peak delays (Fig. 2b,f,i)^^. The observation 
that the slow-speed EPSG is even more sensitive to Klta activity 
is also consistent with our observation that slower EPSPs are 
generally more sensitive to peak shifts (Fig. 2c,e,j). Thus, Klta 
activity promotes inhibition-enforced EPSP peak shifts by 
speeding the decay kinetics of the IPSP, and therefore may 
have a more influential role in facilitating E-I interactions during 
development^"*'^^'^^. 

Inhibition tunes coincidence detection. We have shown that 
inhibition can modulate the peak timing of excitation from one 
side, but it was not known how this influences the integration of 
bilateral excitation. We therefore evaluated the effect of a pre- 
cisely timed IPSP on subthreshold PSP summation and supra- 
threshold action potential probability functions (Fig. 4). Because 
excitation from both sides have equal strength (Supplementary 
Fig. 1) and summate linearly at the soma^' , we injected two 
identical EPSGs at various relative times (Atexc) to simulate ITDs 
in vitro (Fig. 4a, see Methods)^^'^^. We then determined whether 
a single IPSP (timed relative to the contralateral EPSP) could bias 
the excitatory timing of best PSP summation (PSP sum) or 
greatest action potential probability (AP prob). For these 
experiments, it matters at what A^exc the largest summation or 
the most action potentials occur and not when they occur, which 
was the case for previous experiments (Figs 1-3). In the absence 
of inhibition, these functions peaked near zero (PSP sum: 3 ± 2 jis, 
n = 5; Fig. 4b; AP prob: 7 ± 3 jis, n = 7; Fig. 4c). 

Contralateral inhibition modulated coincidence detection 
timing in a manner consistent with its effects on the peak timing 

Figure 3 | The activity of potassium channels facilitates peak shifts. 

(a) Exannple nnodel fittings of EPSPs (top) and IPSPs (bottonn) overlaid with 
corresponding experimentally recorded responses (dotted lines). Scale 
bars, SmV, 5 ms (top) and 2mV (bottom). V.^st- -65mV. (b) Half- 
activation voltage (Vqs, left) and the multiplicative scaling factor for the 
activation time constant (q, right) of the Klta values for each cell in the 
population (n = 16 cells), (c) Comparisons between model-fitted (dark 
markers) and recorded (light markers) EPSP and IPSP amplitude (Amp, left) 
and half-width (HW, right), (d) EPSP peak shifts plotted against Atinh for the 
model prediction (black) and experimental values (grey) for an example cell 
fitting, (e) Same as in d, but the prototypical neuron generated from 
average fit parameters (black) compared with the population of recordings 
(grey), (f) Same as in d, but the Klta ^o.s was depolarized from the 
reference value for the example cell in d (-73mV, light grey) to -53 
(grey) and - 38 (black) mV. The activation kinetics of the Klta was also 
slowed by 100-fold to eliminate Klta activity (turquoise), (g) Same as in f, 
but for the slow-speed EPSG. (h-j) Voltage traces of the EPSPs (h,i) and 
IPSP (j) produced by the model of the cell in (d,f,g) with a Klta ^o.s of - ^3 
(light traces) or - 38 mV (dark traces, left). Normalized (Norm) traces are 
shown on the right to illustrate differences in kinetics. Scale bars, 5 mV, 1 ms 
(h,i), and 2 mV, 1 ms (j). 
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Figure 4 | Inhibition tunes coincidence detection, (a) Schematic 
illustrating three of four inputs examined (left) and timing convention 
(right). Contralateral Inh and ipsilateral excitation are timed relative to 
Contra Exc (Atjnh and Atgxa respectively). (b,d,f,h) Voltage traces for an 
example recording of a subthreshold PSP summation experiment (left). 
Traces were generated from the average-speed EPSG and IPSO and are 
colour-coded to reflect excitatory timing conditions in a. Scale bar, 5mV, 
1 ms. Insets show zooms of the peaks. Summation bias is resolved by 
comparing the peak amplitudes of the flanking 100 |is Ipsilateral-leading 
(yellow) and contralateral-leading (cyan) traces (horizontal dotted lines) 
against zero Atexc (pink). Inset scale bar, 50|iV, 0.1ms. Summated PSP 
amplitudes (PSP sum, normalized to the response of a single EPSP) are 
plotted against Atexc (light) and fitted with Gaussian functions (dark) for 
experiments using the average-(j) and slow (k)-speed EPSG (right). Arrows 
indicate the peak of the function fits (best Afexc)- Negative and positive 
Atexc values indicate Ipsi-leading (l-lead) and Contra-leading (C-lead) 
conditions, respectively. Contra Inh was absent from (b) or occurred at 
Afinh values for average and slow-speed EPSGs that (respectively) advanced 
(0.1/0.2 ms, d), did not substantially shift ( - 0.2/ - 0.1 ms, f) or delayed 
( - 0.6/ - 0.4 ms, h) the peak of the Contra EPSP. Different Atjnh values are 
shown for the average and slow EPSGs because these representative effects 
occurred at different relative timing conditions. n = 5 recordings. (c,e,g,i) 
Same as in b,d,f,h, except that suprathreshold action potential probability 
(AP prob) functions were determined. Scale bar, 5 mV, 1 ms; n = 8-12 
recordings. (j,k) Best Atexc ( ± s.d.) plotted against Atjnh for average-(j) and 
slow (k)-speed EPSGs for PSP summation experiments (open markers) and 
action potential probability experiments (filled markers). 



of contralateral excitation. Inhibitory timing conditions that 
advanced the peak of the contralateral EPSP (Atinh = 0.1 ms) 
biased subthreshold summation and action potential probability 
functions towards ipsilateral-leading excitation (PSP sum: 
-96±8jis, n = 5; Fig. 4d; AP prob: - 73 ± 2 jis, n = 7; 
Fig. 4e). This corresponds to an advanced contralateral EPSP 
that now summates best with an ipsilateral EPSP that occurs 
earlier in time. Timing conditions that did not robustly influence 
the peak timing of the EPSP (Atinh = — 0.2 ms) did not 
substantially bias coincidence detection (PSP sum: — 11 ± 2 |is, 
n = 5; Fig. 4f; AP prob: 6 ± 3 |is, n = 8; Fig. 4g), consistent with a 
different study that was performed with synaptic stimulation^^. 
However, this timing condition (Atinh = — 0.2 ms) did narrow the 
half-widths of the coincidence detection functions (PSP sum: 
387 ±7 versus 467±6|is, P<0.001, n = 5; AP prob: 208 ± 6 
versus 221±8|is, P< 0.001, two-way ANOVA, n = 8), also 
consistent with the synaptic stimulation experiments^^. Finally, 
timing conditions that delayed the peak of the contralateral 
EPSP (Atinh = — 0.6 ms) biased coincidence detection towards 
contralateral-leading excitation, albeit to a lower extent (PSP sum: 
41 ± 6 |is, n = 5; Fig. 4h; AP prob: 41 ± 2 jis, n = 8; Fig. 4i). This 
corresponds to a delayed contralateral EPSP that now summates 
best with an ipsilateral EPSP that occurs later in time. Thus, a 
precisely timed contralateral IPSP can modulate coincidence 
detection in a bidirectional manner, but only under timing 
conditions that had not been previously explored^^. 

Because recent in vivo recordings of MSO neurons revealed 
many cases where the EPSP was much slower than our average 
in vitro measurements^, we also evaluated the slow-speed EPSG 
paired with the average-speed IPSG (Fig. 4b-i,k, chartreuse). The 
peaks of both PSP summation and action potential probability 
functions were biased comparably further towards ipsilateral- 
leading (Atinh = 0.2 ms, PSP sum: — 136 ± 2 jis, n = 5; AP prob: 
— 127 ± 3 jis, n — 7) and contralateral-leading (Atinh — ~ 0.4 ms, 
PSP sum: 64 ± 5 jis, n — 5; AP prob: 66 ± 3 jis, n — 12) excitation. 
This is consistent with the observation that slower EPSPs are 
more sensitive to inhibition- enforced peak shifts (Fig. 2c,e). 
Importantly, the similarity between the peaks of subthreshold and 
suprathreshold coincidence detection functions (Fig. 4j,k) 
indicates that the influence of inhibition on subthreshold 
summation accurately translates into action potential generation. 



Synaptic jitter enhances peak shifts. Because MSO neurons can 
be ITD-sensitive at frequencies between ~100 and ~ 1,500 Hz 
(ref. 48) we returned to our single-side E-I interaction paradigm 
(Figs 1-3) and evaluated the efficacy of inhibitory modulation of 
excitatory timing under such naturalistic conditions. Although 
phase locking is relatively stronger at low frequencies, synaptic 
inputs are actually more jittered . We therefore challenged our 
inhibition- enforced EPSP peak shift paradigm with increasing 
amounts of jitter (Fig. 5). We introduced jitter to four EPSGs and 
IPSGs with gamma distribution functions of up to 640 jis in half- 
width (Fig. 5a,b). The resulting composite EPSPs and IPSPs 
revealed a broadening of event kinetics with increasing jitter, but 
much more pronounced for EPSPs than IPSPs (Fig. 5c). 
Increasing excitatory jitter alone robustly increased peak shifts 
in both directions (for 320 |is excitatory jitter, Atinh = 0.1ms: 
117±8 versus 67±6iis, P<0.001; Atinh =- 0.6 ms: 44±4 
versus 6 ± 3 jis, P = 0.002, two-way ANOVA, n = 6; Fig. 5d,e, 
blue), and increasing inhibitory jitter alone decreased peak shifts 
(for 320 |is inhibitory jitter, Atinh = 0. 1 ms: - 48 ± 4 jis, P = 0.023; 
Atinh = - 0.6 ms: - 5 ± 3 |is, P = 0.344, two-way ANOVA, n = 6; 
Fig. 5d,e, red). Interestingly, increasing both excitatory and 
inhibitory jitter resulted in a net increase in peak shifts (for 320 jis 
excitatory plus inhibitory jitter, Atinh = 0.1ms: —95 ±9 jis. 
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Figure 5 | Synaptic jitter enhances peak shifts, (a) Raster plots showing 
100 trials of four event onset times jittered with gamma distributions of 
increasing function half-widths (left to right). Orange line indicates the 
reference time point (no jitter). Scale bar, 1 ms. (b) Cumulative histogram of 
the onset timing s.d. for the four jittered events, (c) Ten overlaid voltage 
traces for an example recording of the resulting composite EPSPs (top) and 
IPSPs (bottom) generated by jitter functions of increasing half-widths (left 
to right). Kest: -64mV. Scale bar, 2mV, 1 ms. (d) Voltage traces of 25 
trials from the recording in c (from left to right) without jitter, or the 320 [is 
half-width jitter function applied to excitation (Exc) only, inhibition (Inh) 
only, or Exc and Inh. Inhibitory timing conditions that advanced 
(At|nh = 0.1 ms, top) and delayed (Atinh= -0.6 ms, bottom) EPSP peak 
timing are shown. Grey traces show EPSPs and IPSPs alone, and dark traces 
show composite PSPs, aligned in time to the peak of the EPSP alone for each 
trial. Dotted vertical lines and arrow heads indicate the average peak shift of 
all 100 trials in the recording, (e) Average (±s.e.m.) EPSP peak shifts 
plotted against Atjnh for each jitter function applied to Exc only (left), Inh 
only (centre) and Exc plus Inh (right). Without jitter, EPSP peak shifts were 
similar to those in Fig. 1h. Compared with no jitter conditions, introducing 
EPSP jitter (320 |is jitter function half-width) increased peak shifts to 
-200% (Atinh = 0.1 ms, P<0.001; Atinh= -0.6ms, P = 0.002), and 
introducing IPSP jitter reduced peak shifts to ~75% (Atjnh = 0.1 ms, 
P = 0.023; Atinh= -0.6ms, P = 0.344). Introducing EPSP plus IPSP 
jitter increased peak shifts to -150% (At|nh = 0.1 ms, P = 0.014; 
Atjnh = -0.6 ms, P = 0.047). Two-way ANOVA, n = 6 recordings. 



P = 0.014; Atinh = - 0.6 ms: 44 ± 4 jis, P = 0.047, two-way 
ANOVA, n = 6; Fig. 5d,e, magenta). Evaluation of individual 
trials revealed that substantially jittered conditions increase the 
diversity of peak shift magnitudes (Supplementary Fig. 5). 
However, the distribution of all trials (Supplementary Fig. 6) 
shows a trend towards larger peak shifts with equal amounts of 
excitatory and inhibitory jitter. The ability of excitatory jitter to 
outcompete the reduction in peak shifts generated by inhibitory 
jitter can be explained by larger distortion of the EPSP profile 
compared to the IPSP profile with equivalent amounts of jitter 
(Fig. 5c), from which our analysis of synaptic kinetics (Fig. 2) 
predict larger peak shifts. These findings indicate that naturally 
occurring synaptic jitter enhances inhibition-enforced EPSP peak 
shifts, despite the degradation of peak timing precision that was 
recently suggested to diminish its efficacy^ . As ITD tuning at 
low frequencies has been particularly difficult to explain by axon 
length disparity^ ^ these findings support an especially important 
role of inhibition for ITD tuning in the low-frequency range"*^. 



Precise inhibition maintains efficacy at high frequencies. At 

high frequencies, synaptic jitter is much less, but IPSPs naturally 
summate^^. We thus challenged inhibition-enforced EPSP peak 
shifts with 16 pulse trains at 333, 500 and 800 Hz under timing 
conditions that produced a large peak advance (Atinh = 0.1 ms) 
and delay (Atinh = — 0.6 ms) (Fig. 6). Despite substantial 
summation (Fig. 6a, red, and Supplementary Fig. 7), inhibition 
sustained the ability to advance (Fig. 6b, left) and delay (Fig. 6b, 
right) EPSP peak timing, even at 800 Hz (for the 16th event, 
Atinh = 0. 1 ms: - 48 ± 2 jis; Atinh - 0.6 ms: 22 ± 6 |is, n = 8; 
Fig. 6b, bottom). Small changes in peak shifts during the train 
(16th event compared with the first, A^inh = 0.1 ms: — 70 ± 4 |is, 
P<0.001; Atinh = - 0.6 ms: 18±4jis, P = 0.996, two-way 
ANOVA, n = 8; Fig. 6c) can be attributed to frequency- 
dependent interactions between individual IPSPs during 
the train, which alter their effective relative peak timing 
(Supplementary Fig. 7). Inhibition- enforced subthreshold PSP 
summation bias was also similarly maintained at 800 Hz 
(Supplementary Fig. 8), strengthening the translation from 
inhibition- enforced EPSP peak shifts to coincidence detection 
timing. Importantly, at higher frequencies IPSPs cannot be 
discerned during the train regardless of whether they lead or lag 
EPSPs (Fig. 6a, bottom and Supplementary Fig. 8a, bottom). This 
is remarkably consistent with the inability to detect IPSPs from 
in vivo whole- cell recordings^. Thus, inhibition maintains its 
ability to influence EPSP peak timing and synaptic integration at 
high frequencies, despite undergoing substantial summation that 
was recently suggested to diminish its efficacy^^. 



Bilateral inhibition extends the range of tuning. To this point, 
we have considered only the influence of one (contralateral) 
inhibitory pathway on synaptic integration and coincidence 
detection. However, in the full MSO circuit both ipsilateral and 
contralateral inhibition are present^ ^'"^^ We therefore examined 
the additional effect of a second, precisely timed (putative 
ipsilateral) inhibition in biasing subthreshold PSP summation 
(Fig. 7a). In these experiments, we timed the ipsilateral IPSP 
relative to its corresponding EPSP (A^inhi) and used the slow- 
speed EPSG paired with the average-speed IPSG to evaluate the 
extent of inhibitory modulation on synaptic integration under our 
experimental conditions (Fig. 7b). Importantly, an ipsilateral IPSP 
will also interact with the contralateral EPSP, but because it is 
locked in time to the ipsilateral EPSP, the integration of all four 
synaptic inputs can be conceptually simplified as two composite 
PSPs from each side that summate linearly at the soma. 
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Figure 6 | Precise inhibition maintains efficacy at high frequencies, (a) Voltage traces of 16 pulse trains for an example recording of EPSPs (blue), IPSPs 
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individual events that are further analysed in b. Scale bar, 5 mV, 5 ms. The entire train is shown on the left, and a time zoom of the first two events is 
shown on the right. ISI, inter-stimulus interval, (b) The first, second, fourth and sixteenth event for each train in a overlaid for the peak advance (left) and 
delay (right) conditions, aligned in time to the peak of the EPSP alone. Traces are colour-coded to indicate their position in the train (colour code is indicated 
in a, bottom right). IPSPs are omitted for clarity. Insets above each plot are zooms of the peaks, aligned in amplitude. Inset scale bar, 50|iV, 50|is. 
(c) Average (±s.e.m.) inhibition-enforced EPSP peak shifts plotted against each event in the train for peak advance (Atjnh = 0.1 ms) and delay 
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We started with a pre-existing contralateral inhibition that was 
timed (now termed Atinhc) to bias ipsilateral-leading excitation 
(Atinhc = 0-2 ms: — 163±2jis, n — 7; Fig. 7c-e), produce no 
bias (Atinhc= ~ 0.1 ms: — 2±4jis, n = 7; Fig. 7f-h) or bias 
contralateral-leading excitation (Atinhc = — 0-4 ms: 71±7jis, 
n = 7; Fig. 7i-k). In each case, the additional ipsilateral IPSP 
biased the timing of best PSP summation in a manner nearly 
equal and opposite to that of the contralateral IPSP at each Atinhc 
tested. This is consistent with an ipsilateral IPSP that shifts the 
peak timing of the ipsilateral EPSP in the same manner as for the 
contralateral side, which would in fact produce the inverse 
influence on coincidence detection timing (Fig. 8a). Now 
considering the temporal integration of all four synaptic inputs, 
inhibition can bias coincidence detection up to 220 ± 10 jis 
towards ipsilateral-leading excitation (Atinhc = 0-2 and Ati^hi — 
— 0.4 ms; Fig. 7e) and up to 217± 19jis towards contralateral- 
leading excitation (Atinhc = —0.4 and Atinhi = 0.2 ms; Fig. 7k). 
Remarkably, this > 400 |is dynamic range is achieved within a 
timing window of just 0.6 ms (Fig. 8a). This is important, because 
even though we do not know what other internal delays may be 
generated in the circuit, inhibition can provide additional timing 
shifts that largely span the distribution of gerbil ITD functions 
recorded in vivo (Fig. 8b)^~^''*^, which has been thus far 
unsuccessful by axon length disparity-based models of ITD 
tuning^ ^ Interestingly, at timing conditions where inhibition is 
most capable of influencing coincidence detection, IPSPs would 
largely be masked by EPSPs, which is compatible with the in vivo 
observation that even at low frequencies IPSPs remained 
undetected between cycles^. Thus, it is tempting to speculate 
that the strong afterhyperpolarizations following EPSPs recorded 
in vivo^ reflect the re- depolarization phase of a precisely timed 
IPSP that occurs nearly concomitant. 



Discussion 

We demonstrated that precise inhibition modulates coincidence 
detection timing in MSO neurons, supporting a crucial role for 
inhibition in ITD tuning at the cellular level. The cellular basis for 
the near microsecond discrimination of ITDs, an essential 
component of the ability to localize sounds in mammals is 
elusive. As it has recently become clear that axon length disparity 
alone cannot predict the preferred ITD of individual MSO 
neurons^ ^ several alternative models were proposed^'^"*'^^. 

The original inhibition-based model was grounded on the 
finding that blocking glycinergic inhibition in vivo shifted 
contralateral-leading ITD functions towards zero^~^. It was 
thought that a leading contralateral inhibition predominantly 
mediates this effect by delaying the effective peak timing of 
contralateral excitation^'^'^^'^ . This mechanism was questioned 
because it required unusually fast and artificial IPSP kinetics 
and a level of temporal precision that would be diminished by 
synaptic jitter and summation^'^^. We therefore based our 
re-evaluation on measured, mature synaptic responses (Fig. 1) 
under challenging temporal conditions (Figs 5 and 6). We found 
that inhibition was indeed capable of modulating the timing of 
excitation, but that a leading contralateral inhibition alone only 
modestly biased coincidence detection timing towards the 
contralateral-leading side (~50jis; Fig. 4h,i). However, adding 
the ipsilateral source of inhibition^ ^'^^ which has received much 
less attention, could bias coincidence detection a further ~ 150 |is 
towards the contralateral-leading side (Fig. 7k). These data 
strongly support the in vivo pharmacological experiments^'^, but 
demonstrate that the precise timing of inhibition from both sides 
is crucial to explain existing data. 

Our rigorous re-evaluation also revealed that inhibitory control 
over excitatory timing is more diverse than originally proposed. 
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For example, even small differences in synaptic kinetics 
(Fig. 2a-f), the driving force on inhibition (Fig. 2g), the balance 
of excitation and inhibition (Fig. 2g,h), and the activation of Klta 
channels (Fig. 3) markedly influenced the efficacy of inhibition - 
enforced EPSP peak shifts. Moreover, our investigation of a wide 



range of timing parameters revealed conditions that biased best 
coincidence detection towards the ipsilateral-leading side 
(Fig. 4d,e) and others that produced no bias at all (Fig. 4f,g). 
As the actual arrival time of inhibition relative to excitation 
in vivo could be adjusted on a single-neuron basis, inhibition- 
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the difference between the combinatorial bias of bilateral Inh under 
inhibitory timing conditions that are marked in a. Although the magnitude 
of intrinsic delays in the absence of Inh is not known, the range of bilateral 
inhibitory control of coincidence detection timing could modulate the 
preferred ITD of a neuron (best ITD) on an individual basis to an extent that 
exceeds 80% of the population distribution from previously measured ITD 
functions in vivo^'^. Inh, Inhibition. 



enforced ITD tuning in vivo may be more diverse than the 
existing data would suggest^'^. This provides a powerful 
mechanism for tuning individual neurons in the MSO, a central 
issue that remains unexplained by anatomical specializations^ ^'"^^ 
Although inhibition is one of several mechanisms that may work 
in concert to tune neurons to their preferred ITDs, it can 
accomplish this simply by adjustments or adaptations of synaptic 
properties'*^''*'*, a more flexible solution compared with axonal 
length adjustments**. Thus, our findings indicate that the original 
inhibitory model should be expanded. 

Recently, one study evaluated the predictions of the original 
inhibitory model^^, specifically whether a leading contralateral 
inhibition biased ITD functions towards the contralateral side. 
Although there was already evidence that the inhibitory 
input pathways are faster compared with their excitatory 
counterpart^^'^^''*^, this was directly evaluated in a thick slice 
preparation by stimulating the auditory nerve. Inhibition indeed 

10 



preceded excitation, but a single IPSP at the average timing 
condition (A^inh = — 0.3 ms) did not tune ITD functions. It was 
concluded that inhibition, while important for coincidence 
detection, is not relevant for ITD tuning^^. Although this study 
used synaptic stimulation instead of conductance-clamp to 
generate ITD functions, our findings are in accord, as a similar 
timing condition (Atinh = — 0.2 ms) also did not substantially 
influence coincidence detection timing (Fig. 4f,g). Importantly, 
the relative arrival time of excitation and inhibition measured in 
the study were performed in PI 5-20 tissue and thus might not 
represent the fully matured stage of the circuit^^. Moreover, the 
measured relative inhibitory timings spanned a range of — 0.5 to 
— O.lms^^ and are likely even more diverse in vivo, owing 
to cochlear delays and natural activation of the cochlear 
nucleus*^'*^''*^ that are lost in a slice preparation. Therefore, we 
deliberately made no assumptions about the actual arrival times 
of excitation and inhibition and indeed found other timing 
conditions that generated pronounced modulation of coincidence 
detection timing. Finally, the ipsilateral source of inhibition, 
which appeared equivalently as strong as the contralateral source, 
was not previously evaluated in modulating coincidence 
detection^^. Our present results indicate that the ipsilateral 
source of inhibition may be in fact crucial to ITD tuning, and thus 
future studies are required to elucidate the role of its substantial 
input to the MSO. 

The role of synaptic inhibition in the MSO remains a 
topic of intense debate^^"^^. A second recent study investigated 
ITD tuning in vivo using whole-cell recordings. A pronounced 
leading inhibition could not be detected^, contradicting the 
in vitro measurements. It was nevertheless reasoned from the 
lack of detection that the temporal precision of inhibition 
was not sufficient to tune ITD functions. In contrast, our 
findings demonstrate that inhibitory modulation of excitatory 
timing is in fact resilient to synaptic jitter (Fig. 5) and summation 
(Fig. 6). More importantly, under timing conditions where 
inhibition is most effective at influencing coincidence detection 
timing (Fig. 8a), IPSPs would be masked by EPSPs on each 
cycle. Even if inhibition could be detected between cycles at 
low frequencies, its interaction with excitation during high- 
frequency activity would mask its presence, regardless of its 
relative timing (Fig. 6a and Supplementary Fig. 8). Unfortunately, 
as currently available in vivo electrophysiological techniques are 
insufficient to resolve excitatory and inhibitory events 
independently^, the exact arrival time of excitation and 
inhibition to MSO neurons remains unknown. The manner in 
which inhibition tunes ITD functions in vivo thus remains an 
open question. 

The scope of inhibitory mechanisms that modulate ITD coding 
extends far beyond the precise modulation of coincidence 
detection timing we report here. For example, the MSO receives 
GABAergic input, which activates presynaptic GABAb receptors 
on MSO inputs'*^''*^. However, in contrast to phasic glycinergic 
inhibition, GABAb signalling enforces a gain control mechanism 
to adapt ITD sensitivity in dynamic auditory environments^^. 
Furthermore, phasic^^' and sustained^* inhibition are also 
prominent in nucleus laminaris (NL), the functionally equivalent 
avian analogue to the mammalian MSO^^. In contrast to the 
MSO, inhibition in NL is predonimantly GABAergic and imposes 
a shunting and depolarizing effect on NL neurons^^'^^. Here 
phasic and tonic inhibition work in concert to maintain ITD 
sensitivity at different frequencies and intensities^^'^^. The 
inhibitory shunt sharpens the coincidence detection time 
window^ in a frequency region-specific manner, but it does 
not influence the peak timing of incoming excitatory signals^^"^^. 
Thus, both mammals and birds have evolved diverse inhibitory 
mechanisms to modulate ITD processing. 
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Methods 

Tissue preparation. All experiments were performed in accordance with protocols 
approved by the Deutsche Tierschutzgesetz. Male and female Mongolian gerbils 
{Meriones unguiculatus) aged P60-P90 or P17 (Supplementary Fig. 4) were 
anaesthetized with isoflurane [2-chloro-2-(difluoromethoxy)-l,l,l-trifluoro- 
ethane] and decapitated. Brains were removed and placed in a warmed ( ~ 35 °C, 
P60-90) or ice-cold (~4°C, PI 7) dissecting solution containing the following 
(in mM): 200 sucrose, 25 glucose, 25 NaCl, 25 NaHCOj, 3 MgCla, 3 m^o-inositol, 
2.5 KCl, 2 Na-pyruvate, 1.25 NaH2P04, 0.4 ascorbic acid and 0.1 CaCl2 (pH 7.4 
when bubbled with 95% O2 and 5% CO2). Horizontal (110 |im thick, P60-90) or 
coronal (250 |im thick, PI 7) slices were made from the brainstem using a VT1200S 
vibratome (Leica). Slices were then incubated for 30-45 min at 35 °C in a perfusion 
saline containing the following (in mM): 125 NaCl, 25 NaHCOs, 25 glucose, 
3 myo-inositol, 2.5 KCl, 2 Na-pyruvate, 2 CaCl2, 1.25 NaH2P04, 1 MgCl2 and 
0.4 ascorbic acid (pH 7.4 when bubbled with 95% O2 and 5% CO2). 

Electrophysiology. All experiments were performed in perfused saline 
(1 mlmin" ^ and bubbled with 95% O2 and 5% CO2), except that CaCl2 was 
reduced to 1.2 mM. Recording temperature was measured near the slice and 
maintained at 35 ± 1 °C by an SF-28 in-line heater (Warner Instruments) and PH-1 
bath chamber heater (Biomedical Instruments). Tissue was visualized under a 
BX50WI upright microscope (Olympus) equipped with infrared Dodt Gradient 
Contrast optics (Luigs & Neumann), a Polychrome V monochromator (Till 
Photonics) and an Imago camera (TILL Photonics). Whole-cell recordings of 
morphologically identified MSO neurons (Fig. la) were obtained with borosilicate 
glass electrodes using an EPClO/2 patch-clamp amplifier (HEKA Elektronik). 

Voltage-clamp experiments. For voltage-clamp experiments (Supplementary 
Fig. 1), the internal solution contained the following (in mM): 130 Cs- gluconate, 
20 TEA-Cl (tetraethylammonium chloride), 15 HEPES [4-(2-hydroxyethyl)-l- 
piperazineethanesulfonic acid], 5 Cs-EGTA (caesium ethylene glycol tetraacetic 
acid), 5 QX-314 [N-(2,6-dimethylphenylcarbamoylmethyl) triethylammonium 
chloride], 5 Na2-phosphocreatine, 3Mg-ATP, 0.3 Na2-GTP, 0.05 Alexa Fluor 
568 and 0.01 ZD 7288 (4-Ethylphenylamino-l,2-dimethyl-6-methylaminopyr- 
imidinium chloride), adjusted to pH 7.25 and 320 m Osm. Electrodes had 2-3 MQ 
tip resistances, and experiments were performed at — 75 mV (after correcting for 
an estimated liquid junction potential for this internal solution of 15 mV). Series 
resistance was 4-7 MQ and was compensated to a residual of 1.5-2 MQ on the 
amplifier. Experiments were terminated if the uncompensated series resistance 
changed by > 10%. 

Synaptic inputs were stimulated with borosilicate glass electrodes (of 3-4 MQ 
resistance) filled with saline and placed in the vicinity (50-150 |im) of the recorded 
neuron. Ipsilateral and contralateral glutamatergic (excitatory) inputs could be 
independently activated because they are well segregated to the lateral and medial 
sides of the MSO, respectively, (Fig. Ib)^^. As glycinergic (inhibitory) inputs are not 
well segregated in the slice, the population of inhibitory conductance waveforms 
likely reflects inputs from both the lateral and medial nucleus of the trapezoid body 
(lateral and medial nucleus of the trapezoid body)^^. Axons were stimulated with 
brief (0.2 ms) 10-50 V bipolar pulses generated by a Model 2100-isolated pulse 
stimulator (A-M Systems). The stimulus voltage was adjusted such that putative 
single fibres were activated, as indicated by all-or-none threshold responses and 
<25% s.d. in amplitude^ Conductance values (Supplementary Fig. If) were 
calculated from the measured reversal potential for EPSCs (5 mV) and the 
calculated Nernst potential of chloride for IPSCs ( — 44mV). All voltage-clamp 
experiments were performed in the presence of 10|aM R-CPP {3-[(R)-2- 
Carboxypiperazin-4-yl] -propyl- 1-phosphonic acid} and 10|iM SR 95531 [6-Imino- 
3-(4-methoxyphenyl)-l(6H)-pyridazinebutanoic acid hydrobromide] to block 
NMDA receptors and GABA^ receptors, respectively. Glutamatergic and 
glycinergic inputs were isolated with 0.5 |iM strychnine and 20 |iM DNQX (6,7- 
Dinitroquinoxaline-2,3-dione), respectively. 

Single-electrode conductance-clamp experiments. For single-electrode con- 
ductance-clamp experiments (Figs 1, 2 and 4-7 and Supplementary Figs 2 and 
5-8), the internal solution contained the following (in mM): 145 K-gluconate, 
15 HEPES, 5 Na2-phosphocreatine, 3Mg-ATP, 0.3 Na2-GTP and 0.05 Alexa Fluor 
568, adjusted to pH 7.25 and 320 m Osm. Electrodes had 3.5-5 MQ tip resistances, 
and series resistance (4-6 MQ) was 100% balanced on the bridge of the amplifier 
(Supplementary Fig. 4a). The liquid junction potential was estimated to be 15 mV 
and compensated online. Selected synaptic conductance waveforms as measured in 
Supplementary Fig. 1 were delivered to an SM-1 conductance injection amplifier 
(Cambridge Conductance), which calculates instantaneous current commands 
[I{t)] by equation (1). 

/(0=G(0[V(0-£rev] (1) 

This calculation was performed independently for excitatory (reversal potential 
(£j-ev) = 5 mV) and inhibitory (£rev= — 90 mV)^^ synaptic conductance waveform 
templates (G(f)), while simultaneously measuring the membrane potential {V{t)). 
Although single inhibitory fibre conductance was roughly twice that of a single 
excitatory fibre (Supplementary Fig. lb,f), previous studies have indicated that 



there are also about only half the number of inhibitory inputs compared with 
excitatory inputs, and thus the effective inhibitory and excitatory synaptic strength 
is balanced^ Therefore, for subthreshold experiments EPSGs and IPSGs were 
injected with a peak conductance of 30 nS, yielding an E-I ratio of 1:1 (except for 
Fig. 2h, where the IPSG was reduced to 15 nS). To minimize recording time and 
maximize the number of parameters that could be tested in a single recording 
(typically < lOmin)^^, EPSG and IPSG template pairs were randomized and 
then injected at 50 Hz. Current test pulses were delivered directly before and 
after each acquisition epoch to ensure that /h activity had not substantially 
changed^^. Experiments were terminated if membrane input resistance changed 
by > 1 MQ. 

EPSP peak shifts (Figs 1-3, 5 and 6 and Supplementary Figs 3 and 5-7) were 
calculated by comparing the peak EPSP time generated by an EPSG plus IPSG 
injection with that of an EPSG injection alone. Subthreshold PSP summation bias 
(Fig. 4b,d,f,h, Fig. 7, and Supplementary Fig. 8) was determined by Gaussian fitting 
of PSP summation functions. PSP summation functions were generated by 
presenting two identical EPSGs at different relative times (Atgxo in 100 |is 
increments). The amplitude of the composite PSP was then normalized to a 
single EPSP, and the relationship between the normalized composite PSP 
amplitude and Afexc (an in vitro representation of an ITD in the absence of any 
other intrinsic delay mechanisms) was fitted with a Gaussian function. The peak of 
the function fit was defined as the best Afexc- Each timing condition was repeated 
six times in each recording, and analysis was performed on the average voltage 
responses. 

For suprathreshold coincidence detection experiments (Fig. 4c,e,g,i), the timing 
convention was the same as for subthreshold PSP summation experiments, but 
A^exc was presented in 20 |is increments. Because action potential conductance 
thresholds in MSO neurons are extremely sensitive to changes in membrane 
potential, baseline current (192 ± 66 pA, n=l3 recordings) was injected to 
maintain the membrane potential at — 65 mV. Then two identical EPSGs were 
injected with a magnitude (tepsg = 0-3 ms: 43 ± 2 nS, n = l3 recordings; 
T^EPSG = 0-5 ms: 32 ± 2 nS) that brought the neuron to threshold (Gt), defined as 
50% action potential probability at 1 Atexc- The protocol was then repeated 10 times 
at 200 pS (~3%) above Gt, where action potential probability peaked near 100%. 
Repeating these experiments with the inclusion of the IPSG (60 nS, Tipsg =1.5 ms) 
required an adjustment of EPSG magnitude to depolarize the neuron to Gt. 
Adjustments for each Ati^h using the average-speed EPSG (tepsg = 0-3 ms) were as 
follows: 0.1 ms: 50 ± 1 nS, - 0.2 ms: 61 ± 2 nS and - 0.6 ms: 48 ± 1 nS. 
Adjustments for each Afinh using the slow-speed EPSG (tepsg = 0-5 ms) were as 
follows: 0.2 ms: 51 ± 2 nS, - 0.1 ms: 56 ± 2 nS and - 0.4 ms: 50 ± 3 nS. Owing to 
the many trials for these experiments, only 1-3 inhibitory timing conditions were 
performed in each recording. The average action potential probability function 
from all 10 trials in each recording were then averaged for each inhibitory timing 
condition (w = 8-12 recordings per condition). For statistical analysis, the data 
were fitted with Gaussian functions as for subthreshold experiments (Fig. 4b,d,f,h). 

For jitter experiments (Fig. 5), four EPSGs and IPSGs (10 nS each) were jittered 
randomly according to the gamma (F) distribution generated by equation (2). 



The function asymmetry variable a was 40 to mimic the shape of sound-evoked 
spiking output of AVCN neurons^'^, and the function width variable b was 20, 
40, 80, 120 and 160 |is to evaluate different magnitudes of jitter. This generated 
function half-widths of 80, 160, 320, 480 and 640 |is. Iterations of the jittered four 
inputs were presented 100 times at each inhibitory timing condition in each 
recording. 

For train experiments (Fig. 6 and Supplementary Figs 7 and 8), individual 
templates were compounded at 200, 333, 500 and 800 Hz for 16 pulses and then 
injected. Inhibition-enforced EPSP peak shifts (Fig. 6) were calculated relative to 
the excitation -only peak time for each pulse in the train. Analysis for subthreshold 
PSP summation bias (Supplementary Fig. 8) was performed as for single-pulse 
experiments (Figs 4 and 7). Although depression of synaptic currents has been 
reported^ recent observations indicate little if any depression of synaptic potentials 
during ongoing activity^^. Therefore, the conductance amplitude was constant 
for each pulse in the train. This approach thus maximized the effect of summation 
at higher frequencies. Each timing condition was repeated six times in each 
recording. 

Fibre stimulation experiments. For combined fibre stimulation and conductance- 
clamp experiments (Fig. lc,f,h and Supplementary Fig. 2d,e), synaptic conductance 
was injected as for single-electrode conductance-clamp experiments, except that 
an extracellular stimulating electrode was placed on the medial fibre tract, and 
single-fibre stimulation was obtained in current- clamp. Then an EPSG was selected 
to match the decay kinetics of the stimulated EPSP, typically Tepsg = 0.2 or 0.3 ms. 
The strength of the EPSG template was then adjusted to produce an EPSP with 
equal amplitude to the synaptically evoked EPSP (Supplementary Fig. 2d). Then 
the IPSG template amplitude was matched to the EPSG template such that the 
conductance-clamp-simulated EPSP and fibre- stimulated EPSP would have an 
equivalent influence of the IPSP (E-I ratio ~ 1:1). Trials were presented at 10 Hz. 
To eliminate the potential effects of short-term plasticity during the protocol^ each 
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trial was directly preceded by a 20 pulse, 10 Hz train such that the amplitude of 
fibre- stimulated EPSPs was securely at steady stated Each timing condition was 
repeated 10 times in each recording. In separate experiments, single-fibre EPSPs 
and IPSPs were also synaptically evoked in current-clamp for a comparison with 
conductance-clamp- simulated PSPs (Supplementary Fig. 2b,c, blue). 

Dual-electrode conductance-clamp experiments. MSG neuron somata were 
patched sequentially with two electrodes^^, but otherwise performed as for single- 
electrode conductance-clamp experiments. To determine the potential influence of 
hardware errors introduced by injecting conductance and measuring voltage 
through the same electrode, synaptic conductance was always injected through the 
first electrode (E-1), but membrane voltage was measured using the same electrode 
or using the second (passive) electrode (E-2) and compared (Supplementary Fig. 3). 
Both directions were evaluated for each recording. Thus, the terms E-1 and E-2 are 
functionally and not physically defined. 

Data acquisition and analysis. Data were acquired at 100 kHz. Voltage-clamp 
data were additionally filtered at 8 kHz using a shallow three-pole Bessel filter in 
the amplifier. Analysis was performed off-line using Igor Pro (Wavemetrics). 
Example conductance-clamp traces represent an average of 6-10 trials except for 
Fig. 4c,e,g,i and Fig. 5c,d, and Supplementary Fig. 5, which show individual raw 
traces. Example voltage-clamp (Supplementary Fig. lb) and current-clamp 
(Supplementary Fig. 2b, blue) traces represent an average of 30 trials. Group 
averages are presented as mean ± s.e.m., except for the results of Gaussian fittings, 
which are presented as peak time (to) + s.d. Sample sizes for each experiment were 
at least five per group, which generated normal distributions (determined by 
Kolmogorov-Smirnov normality tests, P<0.05). Statistical analysis of the data was 
determined using two-way ANOVA in Prism software (GraphPad). Dunnett's 
post hoc tests were performed to determine statistical significance (P<0.05). 

Computational model. A conductance-based point-neuron model was generated 
to explain subthreshold postsynaptic integration (Fig. 3). The dynamics of the 
membrane voltage V was simulated by current-balance equation (3). 

Cm = - (4 + -^Klta + -^syn + -^leak) (3) 

Cm denotes the measured individual cell capacitance (24.6 pF on average). Ohmic 
currents were assumed in equation (4). 

Uv) =g.d;hi{v-E,) (4) 

gx denotes the peak conductance, and were gating functions with 
cooperativities of m and n, respectively, and Ex was the reversal potential. The 
gating functions adhere to first order kinetics(equations 5 and 6). 



da 



UoQ — a 



d^ 



(5,6) 



floo and /7oo are the steady-state activation, and and x^, are the voltage-dependent 
time constants. The leak current /leak was passive (that is, m = w = 0) with a reversal 
potential Eieak = — 90 mV and a conductance Gieak = 0.8 pS |im ~ ^, estimated from 
the measured remaining input resistance (on average 280 MQ) at — 90 mV when 
Klta and /h channels were blocked with dendrotoxin-k and ZD 7288, respectively. 
The Klta current was modelled with a reversal potential of — 105 mV (ref. 37), an 
average half- activation voltage (Vo.s) of ~ 72.9 mV, and an average activation time 
constant scaling factor (^) of 0.77. This resulted in kinetics for the Klta channels in 
equations (7-10) {V in mV, t in ms, m = 4, n=l). 



floo(V) : 

boo{V) : 



1+6 



6e- 



21.5 

- + 24e^ 



+ 0.35 



1-0.27 



-0.27 



170 



(7,8) 



(9, 10) 



5e- 



M + 10.7 



These equations are, apart from the scaling factor q, the same as given by Mathews 
et al?^ with 35.4 mV added to all voltages. 

The kinetics of the Ih current were modelled by equations 11 and 12. 



aoo(y) - 



,0.1536(V+ 73.97) \ 



1 + e^ 



28.17 + 100.9e 



(11,12) 



3(39-32)/l0 

This is similar to the Ih current of Rothman and Manis^^, but about three times 
slower. The cooperativity was set to m = 1, and the reversal potential was set to 
- 50 mV. 

The synaptic conductance templates were the same as those used for 
conductance-clamp experiments in Fig. 2 (tepsg = 0-3, t:ipsg = 1.5 ms, n = 16; 
''^EPSG = 0-5 ms, n = 5). Conductance-clamp-simulated voltage traces were fitted 



using the following four parameters: Vo.s and q of Klta (Fig- 3b), and the peak 
conductances of 4 and Klta currents. These parameters were adjusted to fit the 
cell- specific PSP kinetics, resting potential and input resistance according to 
conductance-clamp measurements. 
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